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, We study the phase transition in a class of fiber bundle models in which the fiber strengths are 

' distributed randomly within a finite interval and global load sharing is assumed. The dynamics is 

04 ' expressed as recursion relations for the redistribution of the applied stress and the evolution of the 

, surviving fraction of fibers. We show that an irreversible phase transition of second-order occurs, 

^L|i from a phase of partial failure to a phase of total failure, when the initial applied stress just exceeds 

' a critical value. The phase transition is characterised by static and dynamic critical properties. We 

calculate exactly the critical value of the initial stress for three models of this kind, each with a 
different distribution of fiber strengths. We derive the exact expressions for the order parameter, 
the susceptibility to changes in the initial applied sress and the critical relaxation of the surviving 
' (— I ' I fraction of fibers for all the three models. The static and dynamic critical exponents obtained from 

. these expressions are found to be universal. 
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PACS numbers: 46.50.-|-a, 62.20.Mk, 64.60.Ht 



I. INTRODUCTION 



Fiber bundle models describe the collective statics and dynamics of failure in a set of fibers with random strengths 
under the application of a stress (force per fiber) P, . A typical model of this kind is shown schematically in Fig. ^ 
fH , These models are constructed for the purpose of explaining the propagation of fractures in a loaded heterogeneous 
Q ■ material and to determine the conditions under which it breaks completely 0, Q . The latter requires the calculation 
O ] of the strength of the bundle from the strengths of its constituent fibers which^byreasonable assumption, are drawn 
'' at random from a chosen probability distribution P|. In some of the models it is assumed that the load is 

' always divided equally among all intact fibers of the bundle (global load sharing) while in other models |^. 13. lol fiol ITU 
^ . it is assumed that when a fiber breaks, the stress it was last bearing gets distributed only among the fibers next to it 
' (local load sharing). The dynamics or propagation of fracture in a fiber bundle has been characterised in two ways: 
On first, by the probability distribution of bursts of different sizes that occur within the bundle as the stress is gradually 
increased till the bundle breaks completely j^, ilS^ .14, .15,1 ; second, by the of lifetime of a fiber bundle with fatigue 
I under an applied stress It was suggested in |16| that the breakdown of a static fiber bundle with global load 

' sharing can be described as a first-order phase transition, because the surviving fraction of fibers has a discontinuity 
_ at the point of breakdown. However the susceptibility to the applied stress was shown to diver ge a t the breakdown 
. point and for this reason it was later suggested that the transition is of second-order [l^ [l^ |23, |21| . 
' In this paper we report on the universality of the phase transition in a class of fiber bundle models. The main feature 
. of these models is a pair of dynamical recursion relations which, similar to the formulation in [lol l20l | , expresses the 
I ' evolution of the fiber bundle under the application of a finite stress. The strengths of the fibers are assigned randomly 
' ^ [ within a finite interval of values which is true for real fiber bundles. We study three models of this kind: in the first 
5 ■ model the fiber strengths are distributed with uniform density; in the second model the fiber strengths are distributed 
^ ' with a linearly increasing density, which means that there are more strong fibers than weak ones; in the third model 
• • ■ the fiber strengths are distributed with a linearly decreasing density so that there are more weak fibers. From the 
expressions of the fixed points of the dynamics we find that there is a critical initial value of the applied stress in each 
of the three models: on exceeding this critical value the fiber bundle undergoes an irreversible transition from a phase 
^ of partial failure to a phase of total failure. When the initial applied stress is less than or equal to the critical value, 
5^ , only a finite fraction of the fibers breaks as the bundle evolves to a state of mechanical equilibrium; this is the phase of 
partial failure. If, on the other hand, the initial applied stress is greater than the critical value, mechanical equilibrium 
is never reached and the entire fiber bundle eventually breaks down; this is the phase of total failure. We define an 
order parameter which shows that the phase transition is of second-order. As the initial applied stress approaches its 
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critical value from below, the order parameter is found to reduce to zero contimiously following a power- law while the 
susceptibility of the surviving fraction of fibers to changes in the initial stress is found to diverge, also by a power-law. 
We derive asymptotic solutions of the dynamical recursion relations for the surviving fraction of fibers at the critical 
values of the initial applied stress of each model. These solutions show that the critical relaxation of the fiber bundle 
toward the fixed point is a power-law decay. The critical exponents in these power-laws are found to be universal, 
i.e., independent of the distribution of fiber strengths in the bundle. 
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FIG. 1: Schematic diagram of a bundle of clastic fibers with random strengtlis, attached in parallel to a fixed plate at the 
top and a movable plate at the bottom. The fibers are shown as vertical solid lines and the position of the plates by bold 
horizontal lines. The figure shows the fibers arranged in increasing order of strength from left to right; (Tl and an are the 

strengths of the weakest (extreme left) and the strongest (extreme right) fibers respectively, (a) For an initial applied stress 
o"o, o"o < CT-L, the fibers only get stretched from their relaxed position (shown by the bold broken line); the strain 5 measures 
the elastic deformation of the bundle, (b) For aL < ffo < o-r, the fibers begin to break causing a plastic deformation of the 
bundle (now given by the strain S'); the figure schematically shows the strained positions of the broken fibers at the time of 
their breaking in order to give an impression of successive failures. When cro > ctr, all fibers in the bundle break at once. 



II. DYNAMICS OF THE FIBER BUNDLE 



We consider fiber bundle models with global load sharing approximation, i.e., any force applied on a bundle is 

shared equally by all the intact fibers in it. The strength of each fiber is determined by a threshold value (Tthrosh of the 
stress that it can bear, beyond which the fiber breaks. The threshold stress of the fibers in the bundle are distributed 
randomly with a normalised density p(crthresh) within a finite interval [aL,crii] where aL and an are respectively the 
strengths of the weakest and the strongest fiber in the bundle: 



I P(cthresh) dcTthresh = 1- 

O-L 



(1) 
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The probability distribution of the threshold stress is given by: 



{0, < O-thrcsh < CTL 

jj^hrosh p(o-)(lcr, 0-L < CTthrcsh < CTR (2) 
1, O-fl < (Tthrcsh- 

We study the breaking dynamics of the model under the application of an initial stress ctq (for example, by attaching 
an external load to the bottom plate in Fig. ^ where the stress is defined as the force exerted per fiber. The fibers 
whose strengths are less than (Tq break immediately. Following this initial rupture the applied stress gets redistributed 
among the surviving fibers. Since the number of fibers supporting the load has decreased the redistributed stress is 
greater than the initial stress and this causes further breaking of fibers. The process continues till a state of mechanical 
equilibrium is reached, where the surviving fibers have strengths greater than the redistributed stress, or till all fibers 
in the bundle are broken. The state of mechanical equilibrium, if it exists, appears as a fixed point of the model under 
the assigned dynamics. 

The nature of the breaking dynamics allows it to be represented as a recursion relation operating in discrete time- 
steps ^^lEillOl- If Ut is the fraction of fibers in the initial bundle that survive after time-step t, the redistributed 
stress due to global load sharing after t time-steps is 

a. = ^. (3) 
After t + I time-steps the surviving fraction of fibers becomes 



Ut+i = 1 - P{at). (4) 

In a real bundle comprising of a finite number of fibers there will be fluctuations in the local density of fiber strengths 
as well as the load sharing. Such fiuctuatins are ignored in the contruction of Eq. Q and Eq. Q). It follows from 
Eq. (PI and Eq. 10} that the quantities at and Ut evolve by the recurrences 



1 - Pi^t) 



(5) 



and 



Ut+i = 1 - Piao/Ut), Uo - 1, (6) 

which formally define the dynamics of this class of models. With the probability distribution Eq. Q, the recursion 
relation © clearly shows that none of the fibers break for an initial stress ctq < ctl, while on application of (Tq > cn, 
all fibers in the bundle break simultaneously. 

The fixed points of the model, a* and U* , are determined by the relations: 



(7) 



and 



U* + P{ao/U*) = l. (8) 

Though the quantities at and Ut evolve in time till they reach their fixed point values, Eq. (j2Jl shows that their 
product atUt is a constant of motion, always equal to the initial value ctq. All static and dynamic properties of the 
models are consequences of this invariance. 



III. CRITICAL PROPERTIES FOR UNIFORM DENSITY OF FIBER STRENGTHS 



We consider first the case where the random strengths (Xthrcsh of the fibers are distributed with uniform density in 
the interval [aL^aji]. The normalised density function is: 
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0, < (Tthrcsh < CTl 

P(o-thrcsh) = { aa-CTL ' - '^thresh < CT (9) 

0, <7r < CTthrcsh 



and the probability distribution, by the definition in Eq. Q, is given by (Fig. 

0, < CTthrcsh < ctl 



P{<Jthrosh) = { ^"^^^^ < ^thresh < fTfl (10) 

1, ffl < fth rcsh ■ 



We consider only ct^j > ctl as the case of ct^j = ct^ is trivial. 




FIG. 2: (a) The density function p and (b) the probabihty distribution P of random fiber strengths (Jthrcsh distributed with 
uniform density in the interval [(JL,cyR]- In the particular instance shown in the figure ul ~ 0.15 and an — 0.75. 

With this particular choice of the distribution of fiber strengths the recursion relations lO and for an initial 
stress ctl < ctq < ct^j, appear as: 

CTt+l = CTo ( — — — ] (11) 



and 



V CTiJ — CTt 



Ut+i = — - — [<yR-^]. Uo = i. (12) 



The fixed point equations, obtained from Eq. Q and Eq. (ISJ, appear in quadratic form: 



and 



{a*f -ana* +ao{aR-aL) = (13) 



(t/*)2 _ f/* + = 0. (14) 
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Consequently each of the recurrences, Hll() and (|12|l . have two fixed points: 



'1,2 



± {an - ctl) 



1/2 



4((jfl - ctl) 



1/2 



(15) 



TT* 



± 



-,1/2 



CTO 



(16) 



The subscripts 1 and 2 stand for the expressions containing the plus and the minus sign respectively. While 

C/j* are stable fixed points, a\ and are unstable (Fig.OJ. It is clear that the fixed points for the redistributed stress 

and the surviving fraction of fibers are related by: 



U* = 



2,1 



(17) 



(a) 



(b) 



0.75 



D 0-5 




0.25 



0.3 0.375 




0.3 0.375 



FIG. 3: The fixed points of (a) the redistributed stress, and (b) the surviving fraction of fibers for the particular probabihty distribution 
of fiber strengths shown in Fig. |^ In each part of the figure the curve for the stable fixed points is shown by a bold solid line and that for 
the unstable fixed points are shown by a bold broken line. As in Fig. |2]we have = 0.15 and crji = 0.75, so that ctq'''* = 0.234375; the 
position of the critical point is marked by an arrowhead. For ctq < cr^ the fixed points are trivial: since there are no broken fibers cr* = ctq 
and U" =Uq = 1. 



The quadratic equations H13|) and (|14|l show that the initial applied stress ctq has a critical value: 



4(cr_R - ah) 



(18) 



at which their discriminants become zero and only a stable fixed point exists. Since both a and U have physical 
meanings, a state of mechanical equilibrium exists if the quantities (T2 ^''^'^ ^1 positive real- valued; this happens 
only when ctq < cr"'*. For ctq > cr"'*, there are no real- valued fixed points - the dynamics will continue till all fibers 
in the bundle are broken. Therefore a transition ocurs just above cr"'* from a phase of partial failure of the bundle 
in equilibrium to a phase of total failure. The order parameter O for this phase transition is defined in terms of the 
stable fixed point for the surviving fraction of fibers: 
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crit ' 



where 



1 — crit 



(19) 



(20) 



is the value of [/j* under the critical initial stress. As Eq. shows, the order parameter goes to zero continuously 
by a power law as ctq approaches its critical value from below: 



O 







CTO 



1/2 



It is obvious that the critical value must have the lower bound: 



(21) 



^crit 



"0 ^ '^ii 

which, for the expression given in Eq. p8|l . requires that 



(22) 



The above condition (|23() in turn imposes an upper bound on the critical value of the initial stress: 



(23) 



crit ^ CTfl 
^0 



(24) 



Static critical behaviour of the fiber bundle is observed in the susceptibility of the fixed point of the surviving 
fraction of fibers to changes in the initial stress (Tq- From Eq. we see that the susceptibility diverges in the form 
of a power law as the initial applied stress approaches its critical value from below: 



X 



dao 



ctq 



-1/2 



< Co < (To 



(25) 



Dynamical critical behaviour is observed in the process of relaxation of the fiber bundle to a fixed point. At the 
critical point cro = the evolution of the surviving fraction of fibers by II12I) is reduced to the following recurrence: 



Ut+i = 



1 - 



1 



4:i<TR ~ (tl) Ut 



(7o = l, 



(26) 



and its only fixed point, a stable one, is given by Eq. H2U|I . 
The recurrence (|26|) has a closed-form solution: 



Vt - Ul 



-crit 



^l-critl^ ^l-crit) 
l+(l-f/r„ent)(^-l)' 



(27) 



Thus the asymptotic behaviour of the surviving fraction of fibers is a power law decay - a critically slow relaxation 
to the fixed point: 



Ut - C^l-crit 



U- 



1— crit 



cr - (Tl 



t ^ oo 

1 
t' 



(28) 



A special case of this model, with = and aji = 1, was studied in j23|. The critical properties obtained in |2E 
can now be derived easily from the general results of this section. 
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IV. CRITICAL PROPERTIES FOR LINEARLY INCREASING DENSITY OF FIBER STRENGTHS 



We consider next the case where the random strengths (Jthrcsh of the fibers are distributed with hncarly increasing 
density in the interval [o'LtO'ii], a^j > <Ji^. The normahsed density function and the probabihty distribution of the 
fiber strengths are given by (illustrated in Fig. 



and 



P(o'thrcsh) 



rcsh } 



0, 

2(g-th, 



0, 



0, 



(crn-<Ti,)2 



< CTthrosh < 
— CTL < CTthrcsh < CTB. 



rcsh 



< CTthrcsh < CTl 

, aL < CTthrcsh < era 

f_R < tTthresh- 



(29) 



(30) 
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FIG. 4: (a) The density function p and (b) the probabihty distribution P of random fiber strengths (Jthrcsh distributed with 
linearly increasing density in the interval [aL,a_R]. In the particular instance shown in the figure aL = 0.f5 and an — 0.75. 



Here we introduce the transformed quantities: 



To 



0-t 



(^R - 0-L 



(TR. - OL 



(31) 



For an initial stress (Tl < uo < (T_r (or, F^ < Fq < F^ + 1) along with the distribution of fiber strengths given by 
Eq. H3U|) . the recursion relations (jSJ and ^ appear as: 



1 - (r* - TlY 



and 



C/n = 1. 



(32) 



(33) 
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The fixed point equations, |7J and now assume cubic form: 



(r*)^ - 2Tl iT*f + (r| - 1) r* + To = 



(34) 



where F* = <j* /{(tr — ai,), and 



{u*Y + (r^ - 1) {u*Y - {'2TlTo) u* + tI = 0. 



(35) 



Consequently each of the recurrences (|32|) and H33|l have three fixed points - only one in each case is found to be 
stable. For the redistributed stress the fixed points are: 



2 $ 

-Tl + ^K cos-., 



= -ri-Xcos- + VSi^sin-, 
r* = -r^-i^cos-- VSJCsin-, 



(36) 
(37) 
(38) 



where 



(39) 



and 



cos $ = 



(9 - rj) - 27ro/2 
(3 + ri)^^' 



(40) 



Similarly, for the surviving fraction of fibers the fixed points are: 



1 p2 o 
Ul = ^-^ + 2Jcos-, 



1 



— — — J cos h v3j sm — , 



e 



e 



J cos — — V^J sin — 



(41) 
(42) 
(43) 



where 



3 V (^L 



er^ro 



(44) 



and 



cos 9 



1-ri) (ri-i)' + 9riro -27r2/2 



(ri-i) 



3/2 



(45) 



Of these fixed points and C/j* are stable whereas T'l, Fg and C/^ ; ^-^^ unstable (Fig. [SJ. 

Similar to the case in the previous section the discriminants of the cubic equations H34|l and H35|) become zero at a 
critical value u"'* (or, F"'*) of the initial applied stress: 



T^crit 
^ 



0" 


crit 






- (TL 


2 




27 





(3 + Fi)^/^ 



(46) 




FIG. 5: The fixed points of (a) the redistributed stress, and (b) the surviving firaction of fibers for the probability distribution of 
fiber strengths shown in Fig. |1] In each part of the figure the curve for the stable fixed points is shown by a bold solid line and 
those for the unstable fixed points are shown by bold broken lines. We have ol = 0.15 and or = 0.75, so that cr"'* = 0.3375; 
the position of the critical point is marked by an arrowhead. As in the previous case, for (jq < (Jl the fixed points are trivial: 
since there are no broken fibers a* = ao and U* = Uo = 1. 



and then each of the quantities T and U have one stable and one unstable fixed point. As before the critical point 
has the trivial lower bound: 



at' > <yL. (47) 

The expression of F"'* in Eq. (|46ll shows that it approaches the lower bound as F^ — > oo which happens for finite 
values of (Tl and an when {an — a^) ^ 0. It follows that the upper bound for the critical point is also trivial: 



< ^R, (48) 



which is different from the condition (|24|) for the case of uniform distribution. 
At the critical point we get from Eq. H4U|) and Eq. (|45|) : 



cos$ciit = cosecrit = -1 (49) 

or, 

^-crit = 6crit = TT. (50) 

The stable fixed points F2 and f/j* are positive real- valued when Fo < F"'*; thus the fiber bundle always reaches 
a state of mechanical equilibrium after partial failure under an initial applied stress ao < CTq"*. For ctq > cr"'* (or, 
Fq > F"'*), Fj and U* are no longer real-valued and the entire fiber bundle eventually breaks down. The transition 
from the phase of partial failure to the phase of total failure takes place when ao just exceeds cr§'''* and the order 
parameter for this phase transition is defined as in Eq. H19|l : 



(51) 
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Close to the critical point but below it, we can write, from Eq. H45() and Eq. H50|l . that: 



TT- e 



sin 9 

3\/3r™t(3 + r|)3/4(r; 



'Crit 




^ - 1)2 + 6rir^"t]3/2 

and the expressions for the fixed points in Eq. H41|l and Eq. I|42|l reduce to the forms: 



(52) 



and 



where 



-crit ^ _ 1 



2 + 6rLrg"* ^ ° 
rg"*(3 + r2)3/4 



To) 



1/2 



Tj* _ rr* 



crit 



(53) 



(54) 



(55) 



is the stable fixed point value of the surviving fraction of fibers under the critical initial stress cr§"*. Therefore, 
following the definition of the order parameter in Eq. (|51l) we get from the above equation: 



O = 



r™t(3 + r2)3/4 



(ri - 1) 



ePiE- 



crit 



Er 



,1/2 



pcrit 
^ 



(56) 



On replacing the transformed variable Eq by the original ctq, Eq. H56|l shows that the order parameter goes to zero 
continuously following the same power-law as in Eq. 121|l for the previous case when ctq approaches its critical value 
from below. 

Similarly the susceptibility diverges by the same power-law as in Eq. H25|) on approaching the critical point from 
below: 



dC/i* 



dcro 



(r?" - r„) 



-1/2 



To 



-^crit 



(57) 



The critical dynamics of the fiber bundle is given by the asymptotic closed form solution of the recurrence H33|) for 
Eo = Eg"': 



u. 



l-crit/ 



3(E™*) -2EiE™t;7*_ 



t oo, 



(58) 



where E"'* and U^_^j.^^ are given in Eq. H46I) and Eq. (|55|l respectively. This shows that the asymptotic relaxation of 
the surviving fraction of fibers to its stable fixed point under the critical initial stress has the same inverse of time 
form as found in the case of uniform density of fiber strengths [Eq. (|28() ]. 

V. CRITICAL PROPERTIES FOR LINEARLY DECREASING DENSITY OF FIBER STRENGTHS 



Contrary to the case of the previoius section we now consider a fiber bundle where the random threshold values are 
distributed with a linearly decreasing density in the interval [cfl, <7r\, <^r > fi- Instead of Eq. (|29|l and Eq. H30|l we 
now have the following normalised density function and probability distribution (illustrated in Eig.|B|): 



0, < CTthrosh < 

P(crthrosh) = { < CTthrcsh < CTfl (59) 

0, (Tr < CTthrcsh 
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and 



0, < (Jthrcsh < 
P(fTthrcsh) ={ 1 ( "^;":'^-'- )', <JL < fTthrcsh < <JR (60) 

1, ffl < fthrcsh 




FIG. 6: (a) The density function p and (b) the probabihty distribution P of random fiber strengths (Jthrcsh distributed with 
linearly decreasing density in the interval [ctl,o-r]- Similar to the cases shown in Fig. |5|and Fig. |l]we have ctl ~ 0.15 and 
an = 0.75 in this example also. 

With the transformed quantities defined in Eq. (|31|l the recurrences (O and ® for ctl < ctq < ctj? appear as: 

Tt+i = 3 (61) 



and 



Ut+i =i^l + TL-jj-j , Uo^l. (62) 

The fixed point equations are again cubic: 

{r*f -2(1 + Tl) {T*f + (1 + Tl)^ r* - To = 0, (63) 

{U*f - (1 + TLf {U*f + 2(1 + Tl) TqU* - = (64) 
and they have the following solutions: 

n - ^(l + ri) + 2i^'cosy, (65) 

r* = ^(i + rL)-i^'cos|- + y3i^'siny, (66) 

- ^(l + ri)-if'cosy- VSJC'siny, (67) 



12 



where 



l + T, 



3 ' 



(68) 



cos$' 



27ro 

2(i + rr 



- 1 



(69) 



and 



where 



u* 



^^ + 2J cos-, 
^ „ ^' - J' cos — + %/3 J' sin — 



3 

e' 



3 

e' 



J' cos ^ ViJ' sin — 



(70) 
(71) 
(72) 



J' - -v/(i + rL)4-6(i + rL)ro, 



(73) 



cose' 



(i + Ti)^ (i + rL)^-9ro +27r2/2 



(i + rL)^-6(i + rL)ro 



3/2 



(74) 



Here Fg and C/j* are stable fixed points while the rest are unstable (Fig.O. 

The discriminants of Eq. H63|l and Eq. H64|l show that the critical applied stress in this case, erg'''* ' (or, Eg"' '), is 
given by: 



^crit / 



GR - Gh 27 



4 Q 

(i + ri)3 



(75) 



44 



21 (an- a j:)' 

In order to satisfy the condition erg'''* ' > ctl , it requires from Eq. 176|) that 



(76) 



which imposes an upper bound: 



^crit / ^ f^fl 



(77) 



(78) 



Like before, for Fq < Fq'''* ' the stable fixed points are real-valued, which indicates that only partial failure of the 
fiber bundle takes place before a state of mechanical equilibrium is reached; for Fq > F"'* ' the fixed points are not 
real and a phase of total failure exists. The order parameter O of the transition is given by the definition in Eq. I|51|) . 
For Fq = Fg'^'* ' we get the following properties from Eq. ijEHl), Eq. (jZOjl, Eq. (jZU and Eq. ifT^ : 



t^l-crit - t^2-crit - g (1 + Ei) 



(79) 
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FIG. 7: The fixed points of (a) the redistributed stress, and (b) the surviving fraction of fibers for the probabihty distribution 
of fiber strengths shown in Fig. |S| The curve for the stable fixed points is shown by a bold solid line and those for the unstable 
fixed points are shown by bold broken lines. In this example too we have (jl = 0.15 and an = 0.75; here (Jq"^'* = 0.173611, 
marked by an arrowhead. The critical point is located lower than that in Fig.|K|due to abundance of fibers of lower strengths 
compared to the previous case. 



cose;,it = -i or, e;,it = ^ (80) 

and 

cos = ! or, Cit = 0. (81) 

Comparing Eq. H81|) and Eq. H80|l with Eq. (|50|l we sec that the critical values of 8 and Q' are the same whereas those 
of $ and <&' differ by tt radians. 

Near the critical point, but below it, we get from Eq. H7U|) and Eq. H71|) : 



Ul ^ C/r_„.H + ^ (1 + Tl)'/' {Tt' ' T,y'' (82) 

and 

U*^ ^ - ^ (1 + Tl)"^ {Tf' ' ~ Tof'^ . (83) 

Therefore, by the definition of the order parameter in Eq. H51|) and that of the susceptibility in (|57|l . we get in this 

case O oc (To''* ' — Fq)^^^ and x c< (Fo"* ' — To) To Fq''*— . These power laws have the same exponents as the 
corresponding ones in the previous cases and differ from those only in the critical point and the critical amplitude. 

At the critical point the asymptotic relaxation of the surviving fraction of fibers to its stable fixed point [obtained 
as an asymptotic solution to (|62|l ] is again found to be a power law decay similar to Eq. (|28|l and Eq. (|58|l : 



Ut - ui 



4 t/r_crit 

3 



t 



16 
27 



a+TL'f 



(84) 
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The critical behaviour of the models reported in this and the previous two sections show that the power laws 
found here are independent of the form of the probability distribution P. The three probability distributions studied 
have a common feature: the function a* [1 — P (a*)] has a maximum in the interval ((JL,crj^) which corresponds 
to the critical value of the initial applied stress. All probability distributions having this property are therefore 
expected to lead to the same universality class as the three studied here. If the probability distribution does not 
have this property we may not observe a phase transition at all. For example, consider a fiber bundle model with 
P(crthrcsh) = 1 — l/fthrosh, ^thresh > 1- Here (T* [1 — P (cT*)] = 1 and the evolution of the fiber bundle is given by the 
recursion relation Ut+i = Ut/<Jo which implies that there is no dynamics at all for — 1 and an exponential decay 
to complete failure, Ut = (ctq)"*: for fo > 1- There are no critical phenomena and therefore no phase transition. 
However this general conclusion may not be true for finite-sized bundles |24| . 



VI. DISCUSSION 



In this paper we have studied the critical properties of failure in a class of fiber bundle models under an applied 
stress. The models are simple dynamical systems that show an irreversible phase transition. We have determined 
the static and dynamic critical properties associated with the phase transition. Since the models have been defined 
without any fiuctuations in the local density of fiber strengths or in the load sharing, these are equivalent to a mean- 
field theory. We have defined a new order parameter in Eq. H19|l which shows that the transition is of second-order. 
It is supported by facts which arc characteristic of second-order transitions: the susceptibility diverges at the critical 
point and the decay of surviving fraction of fibers with time at the critical point follows a power-law. Besides, the 
ralaxation time of the bundle diverges from both sides of the critical point: for an initial stress infinitesimally below 
the critical point the surviving fraction of fibers will take an infinite time to reach its stable fixed point, whereas 
infinitesimally above the critical point it will take an infinite time to get past the fixed point. By obtaining the same 
critical exponents for three different probability distrbutions of fiber strengths, we are inclined to conclude that these 
exponents are universal. The universality of the critical exponent for susceptibility was also reported in Rcf. 19, 20]. 

The two models studied in sections 4 and 5 seem to be related by the fact that their density functions, Eq. H29|) 
and Eq. (|59|) . can be transformed from one to the other by a reflection on the line (ithrosh = (cl + o'/j)/2 (compare 
Fig. El^a) and Fig. Efa)). But the fixed point equations and their solutions do not have this symmetry. This is 
because the density function p(crthrcsh) does not appear directly in the recursion relations for the dynamics. It is 
the distribution function P((Tthrosh) which appears in the recursion relations. Eq. H3U|I and Eq. (|6Uf) show that the 
distribution functions of these two models are not mutually symmetric about any value of the threshold stress crthrcsh 
(compare Fig.^b) and Fig. Efb)). However a certain relation exists between the critical values of the applied stress 
for a special case of these two models: if = 0, we get from Eq. gni) and Eq. GSJ that o^'^'^Ictr = ^4/27 and 

CTg"* ' /era = 4/27 respectively; therefore we have CTq"' ' /(^r = (o'o"V^-r)^- 

Finally we compare the fiber bundle model studied in this paper with the mean-field Ising model. Though the 
order parameter exponent (equal to i) of this model is identical to that of the mean-field Ising model the two models 
are not in the same universality class. The susceptibility in these models diverge with critical exponents ^ and 1 
respectively on approaching the critical point. The dynamical critical exponents are not the same either: in this fiber 
bundle model the surviving fraction of fibers under the critical applied stress decays toward its stable fixed point as 
t~^, whereas the magnetization of the mean-field Ising model at the critical temperature decays to zero as 
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